Cluster Structure in Cosmological Simulations I: Correlation to 
Observables, Mass Estimates, and Evolution 



Tesla E. Jeltema 1 ' 2 , Eric J. Hallman 3 ' 4 , Jack 0. Burns 4 , and Patrick M. Motl 5 



OO ' ABSTRACT 

O ' 

o 



(N 



We use Enzo, a hybrid Eulerian AMR/N-body code including non- 
gravitational heating and cooling, to explore the morphology of the X-ray gas 
in clusters of galaxies and its evolution in current generation cosmological simu- 
lations. We employ and compare two observationally motivated structure mea- 
sures: power ratios and centroid shift. Overall, the structure of our simulated 
r£3 • clusters compares remarkably well to low-redshift observations, although some 

differences remain that may point to incomplete gas physics. We find no depen- 
dence on cluster structure in the mass-observable scaling relations, Tx — M and 
Yx~ M, when using the true cluster masses. However, estimates of the total mass 
based on the assumption of hydrostatic equilibrium, as assumed in observational 
psj . studies, are systematically low. We show that the hydrostatic mass bias strongly 

correlates with cluster structure and, more weakly, with cluster mass. When the 
hydrostatic masses are used, the mass-observable scaling relations and gas mass 
fractions depend significantly on cluster morphology, and the true relations are 
not recovered even if the most relaxed clusters are used. We show that cluster 
structure, via the power ratios, can be used to effectively correct the hydrostatic 
O ! mass estimates and mass-scaling relations, suggesting that we can calibrate for 

this systematic effect in cosmological studies. Similar to observational studies, 
we find that cluster structure, particularly centroid shift, evolves with redshift. 
This evolution is mild but will lead to additional errors at high redshift. Pro- 
jection along the line of sight leads to significant uncertainty in the structure of 
individual clusters: less than 50% of clusters which appear relaxed in projection 
based on our structure measures are truly relaxed. 
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1. INTRODUCTION 

The growth of large-scale structure in the universe is a powerful probe of cosmology. 
In particular, the largest collapsed structures, clusters of galaxies, are often used as a cos- 
mological probe (see Voit 2005 for a recent review). Clusters are approximately self-similar 
objects, and in the standard paradigm (a cold dark matter universe with structure forming 
from Gaussian initial density perturbations), their bulk properties can largely be predicted 
simply from gravity and the composition of the universe (e.g. the fractions of baryons, dark 
matter, and dark energy). However, as with other cosmological probes (e.g. SN la) the 
accuracy of cosmological studies with clusters is currently limited by systematics, and the 
precise constraints now being sought on parameters like dark energy require a more detailed 
understanding of cluster formation and evolution (e.g. Voit 2005). In particular, cosmologi- 
cal constraints from cluster studies are limited by the accuracy to which we can determine 
cluster mass. 

One major source of uncertainty in cluster mass estimates is the cluster's dynamical 
state. Most methods of estimating cluster mass (e.g. from the velocity dispersion of the 
cluster galaxies, X-ray temperature, or gas mass) assume that the cluster and its gas are 
in equilibrium. However, massive clusters formed relatively recently and are often dynami- 
cally young; clusters form and grow through relatively frequent mergers with other clusters 
and groups of galaxies (e.g. Cohn & White 2005 and references therein). The common 
nature of cluster mergers is supported by the high fraction of clusters observed to have 
significant substructure or disturbed, asymmetric morphologies. The fraction of substruc- 
tured/disturbed clusters observed depends on the method and wavelength considered, but 
for X-ray observations at low redshift this fraction is 40-70% (Mohr et al. 1995; Jones & 
Forman 1999; Schuecker et al. 2001; Kolokotronis et al. 2001). Recent X-ray observations 
with Chandra of clusters out to z ~ 1 also show that cluster structure evolves with redshift: 



clusters are more distu rbed at higher redshifts (IJeltema et al. 20051 : iMaughan et al. 2007 



Hashimoto et al. 20071 ). This evolution is expected in hierarchical structure formation, but 
could lead to additional biases in the selection of clusters and their observable properties 
at high redshift. The evolution of cluster structure in simulations is less clear, but higher 
merger rates are found at high redshift (e.g. Cohn & White 2005; Burns et al. 2007) and 
morphological evolution is detected in a couple of recent simulations (Rahman et al. 2006; 
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Kay et al. 2007). 

Cluster structure is also a powerful probe of the physics underlying cluster formation. 
The formation epoch of clusters and thus the frequency of mergers depends on cosmology 
(e.g. Richstone, Loeb, & Turner 1992). In low density universes, clusters form earlier and 
will be on average more relaxed in the present epoch. In a ACDM universe, the growth of 
clusters slows when the universe transitions to dark energy dominated. Cluster structure 
also depends on gas physics. Cooling, for example, increases the central concentration of 
halos and makes them more spherical (Kazantzidis et al. 2004; Rahman et al. 2006). 

In this paper, we study cluster structure in hydrodynamical si mulations based on Enzo : 



a hybrid Eulerian a daptive mesh refinement (AMR)/N-body code ( INorman fc Bryan 1999 



O'Shea et al. 20051 ). allowing us to directly probe the relationship between cluster structure 
and observable cluster properties. The structure of simulated clusters compared to ob- 
servations tests the reality of the models of cosmology and gas physics employed in current 
simulations. Most previous studies have looked at the effects of mergers on clusters using sim- 
ulations of idealized mergers (e.g. Schindler & Mueller 1993; Roettiger, Burns, & Loken 1996; 
Ricker & Sarazin 2001; Ritchie & Thomas 2002; Poole et al. 2007), individual mergers pulled 
from cosmological simulations (e.g. Mathiesen & Evrard 2001; Rowley, Thomas, & Kay 2004; 
Motl et al. 2005), or simulated clusters selected based on a visually disturbed appearance 
(e.g. Kravtsov, Vikhlinin, & Nagai 2006; Nagai, Vikhlinin, & Kravtsov 2007). These studies 
show that mergers lead to significant variation (usually boosts but it depends on the phase 
of the merger) in cluster luminosity, temperature, and central Sunyaev-Zeldovich (SZ) flux. 
In many cases, mergers are found to boost clusters along the scaling relations rather than in- 
creasing the scatter in these relations, but the exact behavior depends on the type and phase 
of the merger. Promisingly, s ome scaling relat i ons, like Y\ — M and Ysz — M, appear to be 



fairly insensitive to mergers (IMotl et al. 20051 ; iKravtsov et al. 20061 ; iPoole et al. 20071 ). 



In this work, we use observationally motivated measures of cluster structure rather than 
the time evolution of individual mergers to allow a quantitative comparison to the observed 
degree of substructure or asymmetry. The use of a quantitative measure of structure also 
allows us to explore whether these measures can be used either as a tool to select clusters 
or as a third parameter to improve cluster mass estimates. The large volume of these 
simulations (256 Mpc on a side) allows us to study a large number of massive clusters in 
a cosmological setting, giving numbers directly comparable to cluster surveys. We choose to 
focus on simulated X-ray observations of clusters. X-ray observations currently provide both 
the best depth and resolution for studies of cluster structure (compared to the low-resolution 
of weak lensing and SZ observations and the sparse distribution of cluster galaxies), and X- 
ray surveys are the most commonly used to constrain cosmology. We utilize two common 
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measures of structure based on the X-ray surface brightness distribution: the power ratios 
(e.g. Buote & Tsai 1995) and centroid shift (e.g. Mohr et al. 1993), which are described in 
§2 along with the details of the simulations. 

Our results are presented in section 3. We first explore the effect of projection along the 
line of sight on observable cluster structure (§3.1). We then look for correlations between 
cluster morphology and overall cluster properties like mass, luminosity, temperature and Yx 
as well as the effects of cluster structure on the scaling relations used to estimate cluster mass 
(§3.2). In particular, we find a strong correlation between cluster structure and inaccuracy 
in hydrostatic mass estimates. Finally, we compare the structure of our simulated clusters 
to low-redshift observations (§3.3), and look for evolution in cluster structure with redshift 
(§3.4). 



2. METHODOLOGY 
2.1. Numerical Simulations 

In this study, we use numerical galaxy clusters simulated using the cosmological hydro/N- 



body adaptive mesh refinement (AMR) code Enzo ( jO'Shea et al. 20051 ) (http:/ /cosmos. ucsd.edu/enzo) 
Enzo evolves both the dark matter and baryonic fluid in the clusters utilizing the piecewise 
parabolic method (PPM) for the hydrodynamics. The initial dark matter power spectrum 
from Eisenstein & Hu (1999) was initialized on the grid at z = 30. An initial low-resolution 
simulation was used to identify clusters in a volume of 256 h^ 1 Mpc on a side using 128 3 
dark matter particles and grid zones. High-resolution simulations were then performed that 
evolved the entire volume but adaptively refined 50 smaller regions around the largest clus- 
ters identified on the low resolution grid. Each of these 50 regions first is refined with two 
nested static grids, with double the spatial resolution of their parent grid (thus 4 times more 
resolved than the root grid). Within the nested grids, the dark matter particles have a mass 
resolution of 9 x lO 9 /i _1 M . Then, within the nested static grids, we evolve the simulation 
with 5 additional levels of adaptive refinement, again with a factor of two increase in spatial 
resolution at each level. These simulated clusters have a peak resolution of 15.6/i _1 kpc on 
the finest simulation grid. Refinement of high density regions is done as described in Motl at 
al. (2004). The criteria for refinement are both baryonic and dark matter overdensity thresh- 
olds. We assume a concordance ACDM cosmological model with the following parameters: 
Q b = 0.026, Q m = 0.3, Q A = 0.7, h = 0.7, and a s = 0.928. 

At each of 18 redshift outputs, clusters with masses greater than 1O 14 M are identified 
in the refined regions, but unless otherwise noted in this work we chose to use a slightly 
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higher mass cut of 2 x 1O 14 M . We consider only redshifts less than or equal to 1.5, and 
our final sample includes 61 clusters at z = 0.04, 25 clusters at z = 1.0, and 5 clusters at 
z = 1.5. 

The simulations used for this analysis were described in Motl et al. (2005) and Hallman 
et al. (2006), and include models for radiative cooling and star formation with thermal 
feedback from supernovae using the Cen & Ostriker (1992) algorithm. Radiative cooling is 
calcula ted from a tabulated c ooling curve derived from a Raymond- Smith plasma emission 



model (IBrickhouse et al. 19951 ) assuming a constant metallicity of 0.3 relative to solar. The 
cooling curve is truncated below a temperature of 10 4 K. Every timestep, we calculate the 
energy radiated from each cell and remove that amount of energy from the cluster gas 



( IMotl et al. 2004 ; iBurns et al. 20071 ). A full description of the simulated cluster catalog will 



be presented in Hallman et al. (2007). 

We have generated synthetic X-ray images for each cluster in order to apply our structure 
measures to the projected X-ray morphology as has been performed in observational analyses. 
Images are created by calculating a Raymond-Smith X-ray emissivity in the 0.3-7 keV band 
at each grid cell in a 6/i _1 Mpc box centered on each numerical cluster. The values of 
the electron density and gas temperature at each zone are used in this calculation. Then 
the emissivities are integrated along a succession of rays covering the image plane at the 
maximum physical resolution of the cluster. We then convert the physical scale to the 
angular scale that the cluster would subtend at the redshift appropriate to its evolutionary 
state in the simulation. We also enforce a maximum angular resolution of 0.5 arcseconds 
in the images, the maximum resolution of Chandra images. Additionally, in order to more 
directly compare the absolute X-ray flux values of simulated clusters to observed clusters, we 
have rescaled the X-ray surface brightness to account for the difference between the value of 
Qt, used in our simulations and that preferred by the WMAP III results. Images are created 
to include at least the virial radius for each cluster. Clusters are considered to be a single 
object when they overlap within their projected virial radii. 

In the following analysis, X-ray luminosities are determined by integrating the flux 
in the X-ray images over the appropriate aperture and are given for the 0.3-7 keV band. 
Temperatures refer to average spectral temperature and exclude a cool core if present. Cool 
cores are excluded using the temperature profile; we remove regions where the temperature 
decreases by > 20% compared to the surrounding region. Gas masses, including those 
used to calculate Y x , are estimated through deprojection of the X-ray images under the 
assumption of spherical symmetry and using a universal temperature profile (Hallman et 
al. 2006). Since we are interested in the effects of substructure, we do not remove clumps 
from our determination of average cluster properties, as is sometimes done by observers. The 
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removal of substructure in observations will necessarily depend on the depth of the data and 
will not be possible in large shallow surveys, and, in addition, structure along the line of 
sight and shock regions are not removed. We seek to calibrate for these effects rather than 
to remove them. With the substructures included the average temperatures of some clusters 
may be different, most likely lower, than they would be with the substructures removed. 



2.2. Power Ratios 

The power ratios have been described in detail elsewhere (Buote & Tsai 1995, 1996; 
Jeltema et al. 2005); here we give a brief summary of this method. Essentially, this method 
entails calculating the multipole moments of the X-ray surface brightness in a circular aper- 
ture centered on the cluster's centroid. The moments, a m and b m given below, are sensitive 
to asymmetries in the surface brightness distribution and are, therefore, sensitive to sub- 
structure. The physical motivation for this method is that it is related to the multipole 
expansion of the two-dimensional gravitational potential. The multipole expansion of the 
two-dimensional gravitational potential is 



/ i \ 00 i 

V(R,(f)) = -2Ga ln( - J - 2G ^ — — (a m cos m<p + b rn sin rruf>) . 



(1) 



and the moments a m and b m are 



a rn (R) = [ S(f)( J R') m cosm0'dV, 

J R'<R 

b m {R) = [ E(x')(R') m smm<f) , d 2 x', 

JR'<R 



'R'<R 

where x' = (R 1 , 0') and E is the surface mass density. In the case of X-ray studies, X-ray 
surface brightness replaces surface mass density in the calculation of the power ratios. X- 
ray surface brightness is proportional to the gas density squared and generally shows the 
same qualitative structure as the projected mass density, allowing a similar quantitative 
classification of clusters. 

The powers are formed by integrating the magnitude of \l/ m , the rath term in the mul- 
tipole expansion of the potential given in equation (1), over a circle of radius R, 



P m (R) = — / ^ m (R, <t>)^ m {R, <t>)d<t>. (2) 
Ignoring factors of 2G, this gives 

P = [a \n(R)] 2 (3) 
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P m = ^Jj^ K + h 2 m ) • (4) 

Rather than using the powers themselves, we divide by Pq to normalize out flux forming the 
power ratios, P m /P . 

For each cluster in the simulations, we calculate P2/P0, P3/P0, and P4/P0 over a cir- 
cular aperture centered on the centroid of cluster emission (with the origin at the centroid, 
Pi vanishes by definition). P2/P0 measures ellipticity combined with radial fall-off in the 
X-ray surface brightness (P 4 /P is similar but probes smaller scales). P3/P0, as an odd 
moment, is sensitive instead to deviations from mirror symmetry and insensitive to ellip- 
ticity. A high P3/P0 is a clear indication of an asymmetric cluster structure, and this 
ratio shows the stro ngest and most unambiguous evolution with redshift in observations 



( IJeltema et al. 20051 ). In this work, we consider three choices of aperture radius: two radii 
of fixed physical size, 0.5 Mpc and 1 Mpc, for ease of comparison to observations and rsoo, 
the radius at which the mean cluster density is 500 times the critical density. 

Table 1 lists the median power ratios of the simulated clusters for each radius as well 
as the range in power ratios spanned by clusters with the central 80% of power ratios. 
These ranges reflect the overall range of power ratios while eliminating the most extreme 
cases. For the aperture radius of 1 Mpc, we only include clusters whose virial radii are 
greater than 1 Mpc, and it should be noted that this leads to a sample with slightly higher 
average masses than is used for the other two radii. The cluster virial radii all exceed 
0.5 Mpc. The average size of r 50 o for clusters in our sample is 0.8 Mpc, but this radius 
ranges from 0.4 Mpc up to 1.6 Mpc depending on cluster mass and redshift. The overall 
range of power ratios displayed b y the simulated clu s ters compares well to the range of 
power ratios observed in clusters Jjeltema et al. 20051 : iBuote fc Tsai 1996h f see section 3.3 



for further discussion). Investigation of Table 1 reveals that the power ratios of typical 
clusters decrease with increasing radius as enclosed substructure becomes less significant on 
the scale of the aperture size. The same trend is not present for the clusters with the highest 
power ratios due to the fact that increasing radius can lead to the inclusion of additional 
structures within the aperture. 

The power ratios of the simulated clusters are plotted in Figure 1 for both the 0.5 Mpc 
(left) and r 500 (right) aperture radii, and the distributions look similar for R = 1 Mpc. Here 
we show the relations between sets of two of the power ratios. Highly significant correlations 
exist between the power ratios, particularly in the P2/P0 - Pt/Pi plane. S imilar correlations 



are present for observed clusters (IJeltema et al. 20051 : IBuote fc Tsai 19961 ). Such correlations 



imply that the level of structure on different scales increases and decreases together. 
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P3/P0 P3/P0 

Fig. 1.— Correlations among the power ratios for aperture radii of 500 kpc (left) and 
r 500 (right). Clusters are divided in to three roughly equal time bins: z < 0.25 (asterisks), 
0.25 < z < 0.65 (diamonds), and 0.65 < z < 1.5 (crosses). 
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2.3. Centroid Shifts 

Another common measure of observed cluster morphology, based on the first moment of 
the X-ray surface brightness, is centroid shifting (e.g. Mohr et al. 1993; O'Hara et al. 2006). 
In this method the displacement between a cluster's core or peak from the cluster center 
on large scales is used to indicate deviation from dynamical equilibrium. Many different 
methods have been employed to measure centroid shifts. In all methods, the centroid is 
measured in apertures of varying size which could be based on radial steps of set physical 
or angular size (e.g. Mohr et al. 1993), contours of set surface brightness (e.g. O'Hara et 
al. 2006), or radial steps of a set fraction of r 500 (e.g. Poole et al. 2006; Maughan et al. 
2007). The centroid shift is then calculated from the variation of the centroid with radius. 
For example, Mohr et al. (1993) calculated the emission-weighted variation of the centroid 
in circular apertures of set angular size, and Poole et al. (2006) use the RMS offset of the 
centroid from the cluster peak in circular apertures with radii between 0.05r 50 o and r 50 o in 
steps of 0.05r 500 . 

We chose to use a method similar to those used by Poole et al. (2006) and Maughan 
et al. (2007), but here we make a couple of small adjustments. First, we choose to center 
our apertures on the centroid calculated within 0.05r 50 o rather than the cluster peak. The 
cluster peak is difficult to define robustly both for simulated clusters which often have small, 
brightly peaked substructures and for observations where the peak is smeared by noise. In 
practice, we find that the centroid within 0.05r 50 o typically picks out the peak of the main 
cluster. In order to pick the right center for complex cluster morphologies, the centroid 
calculation is begun at large radii with the center allowed to vary and iterated to smaller 
radii to find the centroid within 0.05r 50 o. We then find the centroid in circular apertures 
all centered on this point with radii increasing in steps of 0.05r 50 o out to r 500 . The overall 
centroid shift, (w), is then the standard deviation of the distance between the aperture center 
and the centroid within these apertures normalized by r 50 o- 

Poole et al. (2006) and Maughan et al. (2007) also remove the central 30 kpc surround- 
ing the X-ray peak from the centroid calculation in order to remove the effects of a bright 
cool core and to increase the sensitivity to distortions in the cluster distribution. We also 
investigated the effect of removing the core, but instead we remove the region within our 
smallest aperture, i.e. the region within a radius of 0.05r 50 o- For our simulated sample, which 
encompasses a large range in redshift and mass, 30 kpc represents a varying fraction of each 
cluster. It should be noted that both 30 kpc and 0.05rsoo represent a significantly smaller 
region than the extent of typical cool cores (~ 0.15r5oo; Vikhlinin et al. 2005). However, we 
do not wish to completely remove the cluster core as core structure and central concentration 
relate to cluster dynamical state, but we do remove the very strong peak which could prevent 
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the detection of deviations from symmetry. For the majority of the clusters, the centroid 
shifts are very similar regardless of whether the peaks are removed. 

Table 1 lists the median and ranges of the centroid shifts both with and without the core 
region. Figure 2 shows the relationship between P3/P0 and (w) nocore . The comparison of 
these two structure measures shows a lot of scatter, but they are significantly correlated. For 
comparison, in Table 1 we also list centroid shifts calculated within the physical aperture radii 
of 0.5 Mpc and 1.0 Mpc, which may be more appropriate for the comparison to observations. 
Here the calculation of (w) is quite similar, but the aperture radii are varied in steps of 0.05 
Mpc, and the centroid shifts are normalized to 0.5 Mpc and 1.0 Mpc, respectively. 
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Fig. 2. — Correlation between Ps/Pq and centroid shift (with the cluster peak removed), 
(w) n ocore, for an aperture radius of r 50 o- Clusters are divided in to three roughly equal time 
bins: z < 0.25 (asterisks), 0.25 < z < 0.65 (diamonds), and 0.65 < z < 1.5 (crosses). 
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3. RESULTS 

3.1. Effect of Projection 

Clearly, one important consideration in determining cluster structure from observations 
is projection along the line of sight. Projection can have the effect of either hiding structure 
or projecting unrelated structures to lie close to each other in the plane of the sky. The 
latter case is less concerning in that it can be distinguished observationally with sufficient 
spectroscopy, but this may not be possible for large surveys. In this paper, we typically treat 
clusters which overlap within their projected virial radii as one cluster. With the simulations, 
we create X-ray cluster images for three orthogonal projections and investigate the scatter 
in structure that results for individual clusters. For example, Figure 3 compares the power 
ratios for two projections at both R = 500 kpc (left) and R = r 500 (right). The typical 
deviation in the power ratios between projections is less than a factor of two, which is not 
large for the power ratios (notice that the plots are log-log), but the scatter due to projection 
is significant. The overall variation with projection depends on which ratio and which radius 
are chosen, but for ~90% of clusters the power ratios change by a small fraction (less than 
10%) of the total observed ranges in the power ratios presented in Table 1. A similar scatter 
with projection is seen for the centroid shifts, shown in Figure 4 for an aperture radius of 
r 500 . Again the variation of the centroid shifts with projection for 90% of clusters is less than 
10% of the total ranges in Table 1. 
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Fig. 3. — Scatter in the power ratios due to projection for two orthogonal projections with 
aperture radii of (left) R = 500 kpc and (right) R = r 50 o- The solid line shows a one-to-one 
correlation. 
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Fig. 4. — Scatter in centroid shifts due to projection for two orthogonal projections. Centroid 
shifts are calculated within R = r 50 o both with (left) and without (right) the cluster peak. 
The solid line shows a one-to-one correlation. 
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Interesting questions are how frequently does a disturbed cluster, for example one un- 
dergoing a major merger, appear relaxed in projection. And what fraction of clusters we 
observe to be relaxed is actually disturbed? Clusters exhibit a very large range of mor- 
phologies and structure on many scales, so here we will consider only clusters at the two 
extremes selected based on their power ratios. For each projection, "disturbed" clusters 
are selected based on having P2/P0, P3/P0, or P4/P0 above 1.5 x 1CT 5 , 4.5 x 1CT 7 , and 
1.5 x 1CT 7 , respectively. These limits correspond to the power ratios of observed clusters 
undergoing major mergers and select clusters with roughly the highest 20% of a given power 
ratio in each projection (R = 500 kpc or R = r 500 ). "Relaxed" clusters have P2/P0, P3/P0: 
or P4/-P0 below 1.0 x 10" 6 , 1.5 x 10~ 8 , and 6.5 x 10~ 9 , corresponding to observations of 
relaxed, single clusters and clusters with the lowest ~20% of each ratio. Here we need to be 
careful to account for clusters which may appear disturbed in projection while actually not 
overlaping in 3-dimensions. Using the 3D information from the simulations, we find cases 
where a second cluster (including the region covered by its virial radius) is projected to lie 
within an aperture radius of a given cluster. These cases are then removed from the list of 
disturbed clusters for a given projection. For R = 500 kpc, R = r 500 , and R = 1 Mpc we find 
that 2%, 3%, and 4% of clusters have chance projections with another cluster, respectively, 
although not all of these projections lead to disturbed power ratios. 

We first divide clusters, including all redshifts, in to relaxed and disturbed subsamples 
using the information available from all three projections. Clusters are considered to be 
disturbed if they have power ratios above our limits in one or more of the projections (after 
removing chance projections). For an aperture radius of r 500 , approximately 31% of all 
clusters are disturbed. Relaxed clusters must have power ratios below our limits in all three 
projections, and only about 11% of clusters are relaxed. Of the disturbed clusters, only 
about 55% will look disturbed in any one projection. Of the clusters that appear relaxed 
when viewed in a single projection, about 38% are relaxed from all angles while 8% are 
disturbed (i.e. above our defined power ratio limits) in one of the other two projections. 
The true fraction of relaxed clusters may be even lower, because at the moment of core 
passage during a merger the two cores overlap and the cluster can appear relaxed from all 
angles (Jeltema et al. 2008; Rowley et al. 2004; Poole et al. 2007). Table 2 lists the above 
pecentages for the three aperture radii considered. Percentages are averaged over the three 
power ratios and where appropriate the three projections, but the variations among these 
are not large. In general, these numbers are similar for different radii with the possible trend 
that relaxed clusters are better distinguished at larger radii. 

We ran a similar experiment using the centroid shifts. Table 2 lists the percentages for 
the centroid shifts with the peaks removed, but the numbers are similar if the peaks are 
included. We again define disturbed clusters to be those with the highest 20% of centroid 
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shifts ((w) nocore > 0.047) and relaxed clusters to be those with the lowest 20% of centroid 
shifts {(w) nocore < 0.008). It can be seen from Table 2 that the centroid shifts perform 
a bit better than the power ratios at distinguishing disturbed and relaxed clusters in a 
single projection. For example, only about 4% of clusters which appear relaxed in a single 
projection turn out to be disturbed in another projection, a factor of two better than the 
power ratios. 

Using the 3D information from the simulations it is also possible to determine the 
fraction of clusters which overlap with another cluster within a given radius. Here both 
halos must have masses above 1O 14 M and be sufficiently separated to be identified in the 
simulation as individual halos. We find that 33%, 37%, and 40% of clusters overlap another 
cluster within radii of 0.5 Mpc, r 50 o, and 1 Mpc, respectively. Here we define overlap to 
be when the 3D distance of the center of a second cluster is within a virial radius plus the 
aperture radius. This test does not probe exactly the same types of cluster morphologies 
as the selection of disturbed clusters based on high power ratios or centroid shifts, because 
minor mergers may not result in highly disturbed morphologies and in the later stages of 
mergers the clusters will not be selected as separate halos. It is, however, encouraging that 
both give a similar fraction of merging clusters. 



3.2. Correlations with Other Cluster Properties 

3.2.1. Overall Properties 

We now turn our attention to the correlation between cluster morphology and other 
cluster properties. Do clusters of different masses show more or less structure? If cluster 
mergers affect observable cluster properties like luminosity and temperature, do we observe 
a correlation between these properties and the power ratios? Specifically, we look for cor- 
relations between cluster structure and X-ray luminosity, total mass, X-ray temperature, 
and Yx within a radius of r 50 o- Yx, the product of the gas mass and the temperature, is 
proportional to the total thermal energy of the cluster and is the X-ray equivalent of the 
integrated Sunyaev-Zeldovich (SZ) flux. This measure and the total SZ flux have been shown 
in simula tions to be very l ow scatter proxies for cluster mass and to be fairly insensitive to 



mergers flMotl et al. 20051 : iKravtsov et al. 20061 : IPoole et al. 20071 ). 



Not surprisingly there is quite a bit of scatter in the relations between cluster structure 
and these properties; a cluster of a given mass can have a large range of structures. Only 
cluster luminosity shows a significant correlation with the power ratios. This correlation is 
shown in Figure 5 for P3/P0 versus luminosity with dashed lines at the locations of the limits 
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on P3/P0 defined above for relaxed and disturbed clusters. Table 3 gives the probability of 
a correlation between the three power ratios and each cluster property considered. Proba- 
bilities are derived from a Spearman Rank-Order Correlation test (e.g., Press et al. 1992, 
§14.6), and a probability of one indicates no correlation. The correlation with luminosity is 
in the sense that clusters with higher luminosities tend to have lower power ratios. Since we 
do not find a correlation between mass and the power ratios this trend can not simply be ex- 
plained as a trend for more massive clusters to be more relaxed. The trend of structure with 
luminosity is in the opposite direction of what one might expect given that simulations show 
that mergers can significantly boost cluster luminosity (e.g. Ricker & Sarazin 2001; Rowley 
et al. 2004; Poole et al. 2007). However, relaxed, cooling clusters show very bright cores 
and high luminosities. We also find that during a merger, the cluster luminosity peaks at 
the moment of core passage; during this phase the two cores overlap and the cluster appears 
relaxed (see also Jeltema et al. 2008; Rowley et al. 2004; Poole et al. 2007). We expect 
the core passage phase to be short lived, so this case should be relatively rare, and it may 
be possible to distinguish these cases through temperature variations. P2/P0 also correlates 
with Y x , but the slope is greater than zero at less than 3a significance. 
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Fig. 5. — Correlation between P 3 /P and cluster luminosity within r 50 o- Solid line shows the 
best fit relation. Vertical dashed lines indicate the limits defined in section 3.1 for relaxed 
and disturbed clusters. No significant correlation is seen between the power ratios and total 
mass, temperature, and Y X - 
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Fig. 6. — Correlation between (w) nocore and total mass within r 50 o- Solid line shows the 
best fit relation. Vertical dashed lines indicate the limits defined in section 3.1 for relaxed 
and disturbed clusters. Mass shows the strongest correlation with (w) nocore , but the trend is 
mild. 
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We also look for possible correlations with centroid shift. Table 3 reveals that based 
on the rank-order probabilities there are significant correlations between centroid shift and 
all of the cluster properties except for L x , particularly mass and Y X - However, the best 
fit relations are all fairly flat indicating that these trends are weak. For example, the most 
significant trend is with Yx, shown in Figure 6; clusters with large (w) nocore tend to have 
slightly higher Yx, but again the slope is greater than zero at less than 3a significance. It is 
interesting that we do not see a trend of centroid shift (either (w) or (w) nocore ) with L x , the 
one property that correlates with the power ratios, possibly due to a difference in the types 
of structure probed by these methods. For example, the significant radial dependence in 
the power ratios may make them more sensitive to cool cores. However, in the next section 
we will show that for both the power ratios and the centroid shifts, relaxed clusters have 
luminosities which are high for their masses. We note that the trends presented here also 
depend on the aperture radius chosen. For R = 1 Mpc, we do find significant but mild trends 
in mass and Tx with the power ratios. 

Finally, we note that for any choice of structure measure or radius the change (if present) 
in mass, luminosity, temperature, or Yx with cluster structure is small; the average cluster 
properties change by less than a factor of 1.5 over the full range in cluster structure. Over 
the range of mass considered here (2 x 1O 14 M < M vir < 2 x 1O 15 M ), clusters of different 
masses have a similar range of cluster structures. Below we turn our attention to the effect 
of cluster structure on correlations between different cluster properties, in particular the 
mass-scaling relations. 



3.2.2. Scaling Relations 

Given the well known correlations between cluster properties (e.g. luminosity and tem- 
perature, temperature and mass, etc.), we examine where clusters of different structures fall 
on these scaling relations. For example, Figure 7a shows the relationship between X-ray 
luminosity and temperature within R = r 500 . Here relaxed clusters are shown in red and dis- 
turbed clusters are shown in blue; relaxed (disturbed) clusters have P3/P0 below (above) our 
limits in the one projection plotted. The solid line gives the best fit relation for all clusters 
with the slope fixed at the expected self-similar relation and with the proper evolutionary 
correction applied (i.e. log{E^ x L x ) = 2Mog{T x ) + B where E z = (Q m (l + z) 3 + 1 -Q A ) 1/2 ). 
The dashed and dotted lines give the fits to just the relaxed and disturbed clusters, respec- 
tively. It is clear here that relaxed clusters tend to have higher luminosities for a given 
temperature, and, in fact, our fits show a significant offset between the relaxed clusters and 
the overall cluster population. For the relaxed clusters B = 29.239 ±0.014, while for all of the 
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clusters B = 29.165 ± 0.007. We determine the errors on our fits through 10,000 bootstrap 
trials where clusters are randomly resampled and the fits repeated. The disturbed clusters, 
on the other hand, appear to have a similar distribution in Lx — Tx space to the overall 
population, and the fit to these clusters gives B = 29.155 ± 0.017, very similar to the inter- 
cept for the whole cluster sample. These trends become even more apparent if the selection 
of disturbed and relaxed clusters is done using the information from all three projections, as 
in §3.1. We now define relaxed clusters to be clusters that are relaxed in all three projec- 
tions, and disturbed clusters to be clusters which are disturbed in any one projection after 
removing chance projections. Figure 7b shows the Lx — Tx relation with these selections. 
The disturbed clusters now comprise ~40% of the sample and are consistent with the total 
population, while the relaxed clusters show an even larger offset with B = 29.342 ± 0.027. 

Figure 8 shows the relations between L X) T x , Y x and total cluster mass within R = r 500 
for the same selection of relaxed and disturbed clusters. The L x — M relation shows the 
same trend of relaxed clusters having luminosities that are too high for their masses, while 
the disturbed clusters show no significant offset from the overall population. The Tx — M 
and Yx — M relations, however, show no significant dependence on cluster structure. For 
an aperture radius of 500 kpc, there is a significant trend for the relaxed clusters to have 
temperatures which are low for their masses; however, this trend disappears at larger radii 
like r 50 o and 1 Mpc where the contribution of the cluster core to the temperature becomes 
less important. In contrast, Kravtsov et al. (2006) and Nagai et al. (2007b) find an offset 
of unrelaxed clusters toward cooler temperatures in the T x — M relation. There are several 
differences between our study and theirs, including sample size, average cluster mass, star 
formation and feedback formulation, and selection of disturbed clusters, and it is unclear 
which of these leads to the difference in the derived dependence of the Tx — M relation. 
We do, however, confirm their finding that the Yx — M relation is independent of cluster 
structure. In Figure 8, there are two outliers at low mass representing two redshift outputs 
of a cluster undergoing a major merger where the gas is very offset from the dark matter 
which leads to a problem in centering. This cluster is a known outlier and does not represent 
a general problem. 
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Fig. 7. — L x — T x relation for a radius of r 500 . Top: Red circles show relaxed clusters, blue 
asterisks show disturbed clusters, and the rest of the sample is plotted with diamonds. Here 
the selection of these subsamples was done based on the values of P3/P0 in one projection. 
Middle: The same for relaxed (solid red circles), disturbed (asterisks), and typical (dia- 
monds) clusters with the selection of clusters based on their P3/P0 in all three projections. 
Bottom: The same except that the luminosity of the cluster core within R < 0.15r5oo has 
been removed. Also shown are fits of the form log(E~ 1 Lx) = 2.0log(Tx) + B to all clusters 
(solid line), relaxed clusters (dashed line), and disturbed clusters (dotted line). 
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Fig. 8. — Lx — M {top), Tx — M {middle), and Yx — M {bottom) relations within a radius 
of rsoo- Points show relaxed (solid red circles), disturbed (asterisks), and typical (diamonds) 
clusters selected based on P3/P0 in all three projections. Also shown are fits of the form 
log{y) = Alog{x) + B with the slopes fixed at 4/3 for L x — M, 2/3 for T x — M, and 5/3 for 
Yx — M. Lines show the fits to all clusters (solid line), relaxed clusters (dashed line), and 
disturbed clusters (dotted line). 
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Most of these trends can be explained by the presence of cool cores in the relaxed 
clusters, although core passage events would also create similar offsets. As is frequently 
done in observational studies, we tried removing the cluster cores in the calculation of their 
luminosities. This procedure mitigates the effects of both excess luminosity from cool cores 
and luminosity spikes during the core passage phase of mergers. We adopt the cluster peak 
location from the centroid shift determination (i.e. the centroid within 0.05r5oo) and subtract 
off the l uminosity within R < 0.15r,tsnn from t he total luminosity, similar to observational 



studies ( jVikhlinin et al. 20051 ; iMaughan 20071 ). The results for the Lx — Tx relation are 



shown in Figure 7c. The relaxed clusters are now shifted down in luminosity to have better 
agreement with the overall cluster population, but they still have luminosities which are 
significantly high for their temperatures. The same is true of the Lx — M relation. The 
scatter in the Lx — Tx and Lx — M relations is significantly reduced by removing the 
cluster cores but remains larger than the scatter in the Tx — M and Yx — M relations. On 
the observational side, Hashimoto et al. (2006) find a similar offset between "distorted" 
and "non-distorted" clusters in the Lx — Tx relation when the cores are included, while 
Maughan (2007) finds agreement in the Lx — Yx an d Lx — M relations for relaxed and 
unrelaxed clusters if the cores are removed. Our results confirm that the offset between 
relaxed and unrelaxed clusters as well as the overall scatter can be reduced significantly by 
removing the core, but the remaining offset indicates that, at least in simulations, increased 
luminosity in "relaxed" clusters is also present at larger radii (see also Burns et al. 2007). 

The trends with cluster structure described above are also present, although with typi- 
cally lower significance, if we select disturbed and relaxed clusters based on centroid shift. 



3.2.3. Hydrostatic Masses 



Above we presented the correlations of structure with the true, 3D cluster masses, 
but these are not directly observable. Most often cluster masses are determined from the 
brightness and temperature of the X-ray emitting gas under the assumption that the clus- 
ter is in hydrostatic equilibrium. However, our simulations and previous studies show that 
the hydrostatic masses are systema t ically biased towar d underestimates of the true mass 
( jKay et al. 2004 ; iRasia et al. 20061 ; iNagai et al. 2007a! ). In our simulations, the hydro- 
static masses are lower than the true masses by on average 16%, similar to what is found 
by other authors. Several authors have suggested that this deviation of the hydrostatic 
mass from the true mass is due to additional pressure support from turbulent gas motions 
f lRasia et al. 20041 : IRasia et al. 20061 : iKav et al. 20041 : bolag et al. 20051 : INagai et al. 2007ah 
We might then expect a correlation between this kinetic component and cluster dynamical 
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state. In Figure 9, we show the correlations between the deviation of the hydrostatic masses 
from the true masses versus P3/P0 and (w) nocore . Hydrostatic masses are calculated using the 
3D spherically averaged density and temperature profiles. Our results, therefore, represent 
the best possible result one could get even with perfect data; the typical errors we find in 
the hydorstat ic masses are very s imilar to those found by other investigators using simulated 
observations (jNagai et al. 2007al ). 
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Fig. 9.— Correlation of P3/P0 (top) and (w) nocore (bottom) with the deviation of the 
hydrostatic masses from the true masses within r 500 . Also shown are the best-fit relations. 
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We find a very significant correlation of cluster structure with the inaccuracy of the 
hydrostatic masses, particularly for P^/Pq- The hydrostatic masses of disturbed clusters 
are typically worse underestimates of the true masses with probabilities of no correlation 
of 8 x 10~ 23 and 8 x 10~ 6 for P3/P0 and (w) nocore , respectively. The other power ratios 
show similar trends to P^/Pq. Also shown in Figure 9 are the best fit relations of the 
form AM/M = Alog(X) + B based on a least absolute deviation fit (with the significant 
scatter, fits based on the ordinary least squares estimator appeared to be too biased toward 
outliers). For P 3 /P , A = -0.044 and B = -0.49, and for (w) nocore , A = -0.055 and B = 
—0.26. Similarly, Nagai et al. (2007a) find that the hydrostatic masses of visually unrelaxed 
systems are more biased than visually relaxed systems. Here we quantify the correlation 
with cluster structure, and below we show that this correlation can be used effectively to 
correct the hydrostatic masses. We also find a significant (probability of 6 x 10~ 3 ), but milder 
correlation of mass underestimate with cluster mass, shown in Figure 10. The masses of lower 
mass clusters are typically more underestimated with A = 0.050 and B = —0.88. These 
correlations between deviations from hydrostatic equilibrium and both cluster structure and 
mass can be used as diagnostics for the sources of non-thermal pressure support in clusters. 
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Fig. 10.— Correlation of the true cluster masses with the deviation of the hydrostatic 
masses from the true masses within r 50 o- Also shown is the best-fit relation. 
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The results for cluster scaling relations presented above are encouraging, because they 
show tight relations between observables like Tx and Yx and cluster mass and t hat these rela- 



tions are unbiased wit h respect to cluster structure, confirming previous results (IMotl et al. 2005 



Kravtsov et al. 20061 ). However, the above relations (and previous studies) use the true clus- 
ter masses. As the hydrostatic masses are biased with respect to the true masses, observa- 
tionally normalized Tx — M and Yx — M relations will also be biased, even with the best 
data. Here we investigate the mass scaling relations with the hydrostatic masses and their 
dependence on cluster structure. 

The left side of Figure 11 shows the Lx — Mh y dro, Tx — Mhydro, and Yx — Mhydro scaling 
relations with separate fits to all clusters, relaxed clusters, and disturbed clusters, similar to 
what was shown for the true masses in Figure 8. Here the selection of relaxed and disturbed 
clusters was made using only one projection, corresponding to what is actually observable. 
These fits are listed in Table 4 as well as the fits to the total cluster sample using the true 
masses which we will take as the "true" scaling relations. Again, the slopes have been fixed 
at their predicted self-similar values (4/3 for Lx — M, 2/3 for Tx — M, and 5/3 for Yx — M). 

It is clear from Figure 11 that using the hydrostatic masses leads to both a significant 
increase in the scatter of these relations as well as systematic offsets due to cluster structure. 
We find significant offsets in the Tx — Mhydro and Yx — Mhydro relations between relaxed and 
disturbed clusters (the known offset of the relaxed clusters in the Lx — M relation is less 
obvious here). In addition, the mass scaling relations are significantly different than those 
obtained using the true masses even if the most relaxed clusters, which have the least biased 
hydrostatic masses, are used. For example, the normalization of the Yx — Mhydro for the 
relaxed subsample is 27% higher than the relation for all clusters using the true masses (11a 
significance). This offset easily explains the offset found by Kravtsov et al. (2006) between 
simulated and observed clusters (see also Nagai et al. 2007b). This result implies that even 
if deep observations of relaxed clusters are used to normalize the mass scaling relations the 
normalizations will be systematically biased. 



Fig. 11— Left: L x - M hydro (top), T x - M hydro (middle), and Y x - M hydro (bottom) 
relations where the masses are calculated assuming hydrostatic equilibrium and using the 
3D spherically averaged density and temperature profiles. All quantities are calculated within 
a radius of r 50 o- Red circles show relaxed clusters, blue asterisks show disturbed clusters, 
and the rest of the sample is plotted with diamonds. The selection of these subsamples was 
done based on the values of P3/P0 m one projection. Lines show the fits with the slopes 
fixed as in Figure 7 to all clusters (solid line), relaxed clusters (dashed line), and disturbed 
clusters (dotted line). Right: Same except that the hydrostatic masses have been corrected 
based on P3/P0 for each cluster. 
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However, the strong correlation between cluster structure and the deviation of the hy- 
drostatic masses from the true masses shown in Figure 9 means that we can use cluster 
structure to at least partially correct the hydrostatic mass estimates. The right side of Fig- 
ure 11 shows the scaling relations that result if we use the fit in Figure 9a to adjust the 
hydrostatic cluster masses based on their P^/Pq. The corresponding fits are listed in Table 
4. The Tx — M and Yx — M relations are now unbiased with respect to cluster structure 
showing no significant offsets for the relaxed and disturbed subsamples. The significant offset 
of relaxed clusters in the L x — M relation present for the true masses is now also reproduced. 
Even better, the normalizations of all of the scaling relations are now in very good agree- 
ment with the relations based on the true masses, although the scatter about these relations 
remains higher. 

The above results are for an assumed slope; however, we also find that the best fit 
slopes of the mass scaling relations change if we use the hydrostatic masses. Figure 12 and 
Table 5 compare fits with variable slope to the mass scaling relations (for all clusters) using 
the true, hydrostatic, and corrected masses. These fits use the bisector modification to the 
BCES method in Akritas & Bershady (1996), a regression method commonly employed by 
observers because it allows for intrinsic scatter and non-uniform measurement errors in both 
variables. Here it can be seen that the scaling relations based on the hydrostatic masses do 
differ significantly in both slope and normalization from the true relations. Correcting the 
hydrostatic masses based on cluster structure changes the best-fit slopes only slightly, but 
does shift the normalization to make the relations in overall better agreement with most of 
the clusters. The lack of change in slope after mass correction is unsurprising given the lack 
of correlation between mass and cluster morphology shown in section 3.2.1. Correcting for 
the additional correlation between mass and the bias in the hydrostatic mass (Figure 10) 
would bring the slopes in to better agreement. 

Table 5 also lists the scatter about the fit relations in terms of the scatter in log(X) 
(mass) for fixed observable. Employing the hydrostatic masses significantly increases the 
scatter, and the scatter is only slightly reduced by correcting the masses. We have shown 
that we can eliminate systematic biases in the mass-scaling relations; the scatter remains 
higher than the true relations due to the scatter in the relationship between our structure 
measures and deviations from hydrostatic equilibrium. As we showed in section 3.1, one of 
the major sources of scatter in the relationship between observed structure and dynamical 
state is projection, and the nature of this scatter can be better understood through studies 
of this kind. Similar to Kravtsov et al. (2006), we find the lowest scatter in mass for the 
Yx — M relation. 
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Fig. 12. — Left: Comparison of the fits to the T x — M (top) and Y x — M (bottom) relations 
using the true masses (solid line), hydrostatic masses (dashed line), and corrected masses 
(dotted line). For comparison, the points show the properties of the simulated clusters and 
their true masses. Right: The solid line shows the difference between the masses that would 
be estimated from the Tx — Mhydro (top) and Yx — Mh y dro (bottom) relations and the masses 
derived from Tx — M true and Yx — M true . This difference is presented as an absolute fractional 
deviation. Also shown are the absolute fractional deviations relative to the true relations if 
the hydrostatic masses are corrected by P$/Pq (long dashed), by (w) nocore (dot dot dashed), 
and by the median error in the hydrostatic masses (dotted). 
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The advantage of correcting the hydrostatic masses using P3/P0 is shown clearly on 
the right side of Figure 12. The solid line shows the absolute fractional difference between 
the masses that would be estimated from the T x — M hydro and Y x — M hydro relations and 
the masses derived from T x — M true and Y x — M true , while the dashed line shows the same 
difference after the hydrostatic masses are corrected by P^/Pq. The corrected masses are 
closer to the true masses over most of the ranges in T x and Y x , except at the highest 
masses. At most the corrected masses are off by 20% compared to the true masses, and 
unlike the hydrostatic masses, which for the most part show a systematic shift to lower 
masses, these can be either over or underestimates. The same result can be seen if we 
compare the corrected masses for individual clusters: the corrected masses are closer to the 
true masses than the hydrostatic masses for about 75% of the simulated clusters. Also shown 
in Figure 12 are the resulting mass errors if we instead correct the hydrostatic masses based 
on centroid shift (dot dot dashed line) or based simply on the median error in the hydrostatic 
masses (dotted line). The centroid shifts show a less significant correlation with mass error, 
and this correction performs worse than the power ratios over the full ranges of T x and Y x . 
It appears that while the centroid shifts are more robust against projection than the power 
ratios, they are less sensitive to deviations from hydrostatic equilibrium. Correcting simply 
by the median error in the hydrostatic masses of course leads to an improvement in the mass 
estimates but performs worse than or similar to P3/P0 everywhere, showing the need for a 
structure dependent factor. 

Here we have not explicitly shown the effect of cluster structure on gas mass fraction, 
another important cosmological tool. A couple of recent studies have shown that gas mass 
can be es timated fairly accura t ely (within ~ 10%) with a relatively small bias for unrelaxed 



clusters ( jHallman et al. 20061 ; iNagai et al. 2007al ). In this work, we consider gas masses 



estimated using a universal temperature profile and deprojection under the assumption of 
spherical symmetry, and we find a significant correlation between cluster structure and devi- 
ation of the estimated gas masses from the true gas masses. The errors in gas mass, though, 
are fairly small. The gas masses of relaxed clusters are biased high by less than 5% while the 
gas masses of disturbed clusters are high by ~10%. However, the error in the hydrostatic 
masses and its structure dependence will translate directly in to error in the gas mass frac- 
tions. Even considering only our relaxed subsample, the typical error in hydrostatic mass is 
12% and individual "relaxed" clusters have errors as large as 56%. For disturbed clusters, 
the ~10% overestimate in gas mass will add to the larger underestimate in hydrostatic mass. 
A similar result was found by Nagai et al. (2007a) for visually unrelaxed systems. 
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3.3. Comparison to Observations: Low Redshift 

An important question is how well the simulations reproduce the observed distribution 
of cluster morphologies. Cluster structure and its evolution depend on both cosmology and 
gas physics and can act as a sensitive test of how well current simulations are reproducing 
the observed universe. A detailed comparison to observations including the effects of noise, 
instrumental response, and selection is beyond the scope of the present paper but will be 
explored in future publications. Here we limit our scope to a broad comparison of the 
structure of simulated clusters to low-redshift observations where the signal-to-noise is high 
and observational effects can be expected to be minimal. 

We compare the power ratios of our simulated clusters to the ROSAT observed sample 
from Buote & Tsai (1996). We limit both samples to include only clusters with redshifts 
less than or equal to 0.1 and luminosities greater than 10 44 ergs s -1 . For the simulations 
this includes three redshift bins at z — 0.04, 0.08, and 0.1. At these redshifts, the ROSAT 
resolution (~ 20") is reasonable compared to the resolution of the simulations (28" at z=0.04 
and 12" at z=0.1), and the ROSAT field-of-view, unlike Chandra's, is large enough to study 
cluster structure to reasonably large radii. For direct comparison to the observations we 
consider the power ratios within the physical aperture radii of 0.5 Mpc and 1 Mpc. Table 
6 lists the median power ratios for these two radii. We also list the probabilities from a 
Wilcoxon rank-sum test (e.g., Walpole & Myers 1993) that the mean power ratios of the 
observed and simulated samples are the same and the probabilities from a Kolmogorov- 
Smirnov test (e.g., Press et al. 1992, §14.3) that the distributions of power ratios are the 
same. Figure 13 compares the distributions of power ratios in the simulations (solid line) 
and observations (dashed line) for R = 0.5 Mpc. 
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Fig. 13. — Comparison of the distributions of power ratios within R = 0.5 Mpc in simulations 
(solid line) and observations (dashed line) for low-redshift clusters (z < 0.1) with luminosities 
greater than 10 44 ergs s _1 . The number of observed clusters (42) is normalized to give the 
same total number of clusters as in the simulations. The observational sample is taken from 
the ROSAT sample of Buote & Tsai (1996). 
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These comparisons show that clusters in our simulations have remarkably similar struc- 
tures to observed clusters compared to what might have been expected. The power ratio 
distributions have similar peaks, and we find no statistically significant differences in either 
the mean power ratios or the distributions of the two samples (with the possible exception 
of P2/P0 at R = 0.5 Mpc). However, investigation of Figure 13 reveals that the simulations 
show a broader distribution of power ratios than observed, particularly for the even power 
ratios. The simulations appear to include both more clusters with very high power ratios and 
a tail of clusters with low power ratios extending to lower power ratios than are observed. 
Investigation of those clusters in the simulations which have power ratios lower than what 
is observed reveal that these are typically low mass systems with luminosities that are high 
for their masses. Clusters with power ratios higher than observed, on the other hand, tend 
to have larger masses. For example, Figure 14 shows the Lx — M relation for the simulated 
clusters in Figure 13 with very low power ratio (circles) and very high power ratio (squares) 
clusters highlighted. 
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Fig. 14.— L x — M relation (R = 0.5 Mpc) for simulated clusters with z < 0.1 and 
luminosities greater than 10 44 ergs s -1 . Solid circles mark clusters with power ratios which 
are low compared to observed clusters, and solid squares mark clusters with power ratios 
which are high compared to observed clusters. 
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Figure 14 shows a similar trend for low power ratio clusters to be offset to higher 
luminosities as was seen in section 3.2.2, but here we use a stricter cut-off on morphology. 
These observations strongly suggests the role of cool cores in low-mass clusters in producing 
at least the low power ratio tail of clusters. In fact, all but one of the eighteen low power 
ratio clusters marked in Figure 14 are classified as a cool cores by the criteria defined in 
Burns et al. (2007). The radial dependence of the power ratios means that if clusters are 
more strongly peaked in the simulations compared to observations, the power ratios of single 
component clusters will be significantly lower. At the high end, small but bright, strongly 
peaked substructures, bright shocks, or cold fronts could contribute to increasing the power 
ratios. On the other hand, the median P2/P0 and P4/P0 of clusters in the simulations are 
slightly higher than observed. Together with the fact that P3/P0 shows the best match to 
observations, this suggests that the simulated clusters are on average a bit more elliptical than 
observed. As cooling acts to reduce cluster ellipticity (Kazantzidis et al. 2004; Rahman et 
al. 2006), our results could indicate the need for additional cooling in more massive clusters. 
A couple of other studies have also found an offset between the observed and simulated 
ellipticites of clusters, but not in a consistent direction (Rahman et al. 2006; Flores et 
al. 2007). This comparison likely depends on the specifics of the cooling and feedback 
prescriptions used. 

Examples of clusters in the simulations with very low and very high power ratios are 
shown in Figures 15 and 16, respectively. The three clusters with very small power ratios 
in Figure 15 are all low mass, cool core clusters with very strong surface brightness peaks 
and large dips in temperature in the core. The clusters in Figure 16 with high power ratios 
are unsurprisingly strongly merging with one or more bright substructures, shocks in both 
the surface brightness and temperature maps, and no cool core. These clusters are viewed 
near the peak of the merger. The third cluster, in particular, has a main component which 
is very disrupted showing no clear core. 
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Fig. 16. — Examples of clusters in the simulations with very high power ratios, higher than 
typically seen in observations. Images are 2 Mpc on a side. 
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Observational effects which might contribute to the observed cluster morphology include 
resolution and noise, both of which could smooth out bright structures. We tested both of 
these effects using the simulated clusters shown in Figures 15 and 16. First the images were 
smoothed with a Gaussian filter with width equivalent to the ROSAT resolution. However, 
even at z=0.1, the ROSAT resolution is within a factor of two of the resolution of the 
simulations and sufficiently small to resolve the typical scale of structure probed by the 
power ratios, and this smoothing has little effect. We then added noise appropriate for a 10 
ksec observation with ROSAT, which is the typical exposure time of the clusters in the Buote 
& Tsai (1996) sample. The addition of noise has almost no effect on the high power ratio 
clusters in Figure 16, which is not surprising. Jeltema et al. (2005) showed that for clusters 
with significant structure, the power due to noise is a very small fraction of the total power. 
Noise also has little effect on P 2 /Po, because it is the lowest order power ratio; however, this 
is the ratio that shows the largest discrepancy between observations and simulations. P3/P0 
and P4/P0 for the relaxed clusters in Figure 15 do, for the most part, increase enough with 
the addition of noise to be consistent with the values of these ratios for relaxed clusters in 
the observed sample. 

The above tests reveal that observational effects like noise can account for some dif- 
ference in the power ratio distributions between observations and simulations, but that 
"gastrophysics" also plays a role. It should be emphasized that the agreement in cluster 
structure between simulations and observations, at least at low redshift, is remarkably good 
with no statistically significant differences in either the means or distributions of the power 
ratios (except perhaps a mild offset in P 2 /P ). The agreement indicates the general success 
of both the underlying cosmological model and the gas physics employed here. The disagree- 
ments (like the broader distribution of power ratios in the simulations seen in Figure 13), on 
the other hand, suggest that with a more detailed treatment of observational effects cluster 
structure will provide a sensitive test of "gastrophysics" in simulations. For example, our 
results suggest the need for a reduction in the peakiness of low mass cool cores but an in- 
crease in cooling in high mass clusters. We may need an additional feedback mechanism like 
AGN or a slightly different balance between the cooling and heating prescriptions. Another 
general success of these simulations is that they are the first set of simulations to produce 
both cool core and non-cool core clusters in the same volume (Burns et al. 2007). However, 
the fraction of cool cores is low compared to observations of X-ray selected clusters, and 
there are very few high-mass cool cores in our sample (Burns et al. 2007), supporting our 
conclusions based on cluster structure. 

In this section, we do not consider centroid shifts, because no low redshift observed 
sample of clusters using the same formulation of centroid shifts exists. We note that our 
centroid shifts are very similar in range and median to those found by O'Hara et al. (2006), 
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who define centroids within surface brightness contours rather than at at fixed radii. For 45, 
z < 0.2 clusters, they find (w) between roughly 0.005 and 0.1 r 500 with a median value of 
~ 0.02. 



3.4. Evolution of Structure 

Jeltema et al. (2005) showed using Chandra observations that cluster structure evolves 
with redshift: high-redshift clusters are on average more disturbed than low-redshift clusters. 
Specifically, the average P3/P0 power ratio increases with redshift. The evolution of cluster 
structure has now been confirmed by Maughan et al. (2007) using centroid shifts and by 
Hashimoto et al. (2007) using asymmetry in the X-ray surface brightness distribution. A 
direct comparison to these observations, particularly at high redshift, should include a more 
detailed treatment of selection and observational uncertainties, and we will defer an in depth 
discussion of this to a future paper. In this section, we focus on whether or not cluster 
structure is observed to evolve in our simulations. 

To quantify the significance of possible evolution in structure, we divided the clusters 
into three redshift bins of relatively equal time intervals: low redshifts (z < 0.25), inter- 
mediate redshifts (0.25 < z < 0.65), and high redshifts (0.65 < z < 1.5). We then used a 
rank-sum test to find the probability that the average power ratios and centroid shifts in 
these redshift bins are the same. The results are given in Table 7. For the power ratios, we 
list the results for both R = and, for better comparison to observations, R = 0.5 Mpc. 
We find a significant increase in the average centroid shifts with redshift, but no significant 
evolution in the power ratios. For example, Figure 17 shows the median P3/P0 (top) and 
(w) n ocore (bottom) within R = r 50 o for each redshift bin in the simulations. For compari- 
son, the dashed lines indicate the lowest 20% and highest 20% of each structure measure 
versus redshift. The median centroid shift does appear to increase slightly with redshift, 
while P3/P0 shows little or no evolution. The right side of Figure 17 shows histograms of 
the distributions of P3/P0 and (w) nocore in the three redshift bins considered above; here the 
evolution in (w) nocore can also be seen. These plots show, however, that any evolution in 
structure is mild compared to the range of cluster morphologies present at any redshift (i.e. 
there are relaxed and disturbed clusters at all redshifts). 
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Fig. 17.— Left: Median P3/P0 (top) and (w) nocore (bottom) within R = r 50Q for each 
redshift bin in the simulations. Also plotted are dashed lines indicating the lowest 20% 
and highest 20% of each structure measure versus redshift. Right: Histograms showing the 
distribution of P3/P0 (top) and (w) nocore (bottom) in three redshift bins: z < 0.25 (solid 
line), 0.25 < z < 0.65 (dashed), and 0.65 < z < 1.5 (dotted line). 
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Our simulated sample, however, is mass selected, while the observational samples detect- 
ing structure evolution are primarily composed of clusters from flux-limited X-ray surveys. 
Jeltema et al. (2005) used a luminosity limit of Lx > 2 x 10 44 ergs s _1 in their sample selec- 
tion, and for comparison we also list in Table 7 the probability of evolution in the simulations 
for this sample selection. With a luminosity rather than mass cut, we find significant evo- 
lution in all three power ratios and the centroid shifts for R = r 50 o- This evolution appears 
to be mostly due to an increase in the average cluster structure from low to intermediate 
redshifts. The stronger evolution of the power ratios in a luminosity-selected sample may be 
due to a preferential selection of relaxed, cool core systems (§3.2) which show more significant 
cool cores at low redshifts. For an aperture radius of R = 0.5 Mpc, the same radius used by 
Jeltema et al. (2005), the power ratios do not show significant evolution (with the exception 
of P3/P0 between low and intermediate redshifts), but the use of a constant physical radius 
aperture could act to hide a mild evolution in cluster structure due to the fact that cluster 
virial radius decreases with increasing redshift. 

Figure 18 compares the evolutions in P3/P0 (top) and (w) nocore (bottom) in the simula- 
tions to the current observational samples. Here we use a luminosity cut of Lx > 2 x 10 
ergs s _1 for both the observations and the simulations and appropriate choices of aperture 
radii to match the observations (R = 0.5 Mpc for P3/P0 and R = r 50 o for (w) nocore ). The 
observed power ratios are taken from Buote & Tsai (1996) for z < 0.1 (lowest-redshift bin) 
and Jeltema et al. (2005) for higher redshifts (higher-redshift bins). The observed centroid 
shifts are taken from Maughan et al. (2007). The observational data points show the median 
value in redshift bins containing at least 20 clusters, with the exception of the highest-redshift 
P3/P0 bin which has 15 clusters. As stated above, in this paper we will not seek a detailed 
quantitative comparison between structure evolution in simulations and observations, but 
since significant evolution in cluster structure is detected in observations, here we seek a 
qualitative comparison. 
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Fig. 18. — Evolutions in P3/P0 (top) and (w) nocore (bottom) in the simulations compared 
to the current observational samples. We use a luminosity cut of Lx > 2 x 10 44 ergs s _1 for 
both the observations and the simulations and appropriate choices of aperture radii to match 
the observations (R = 0.5 Mpc for P3/P0 and R = rsoo for (w) nocore ). The observed power 
ratios are taken from Buote & Tsai (1996) for z < 0.1 (lowest-redshift bin) and Jeltema et al. 
(2005) for higher redshifts (higher- redshift bins). The observed centroid shifts are taken from 
Maughan et al. (2007). The observational data points (asterisks) show the median value in 
redshift bins containing at least 20 clusters, with the exception of the highest-redshift P3/P0 
bin which has 15 clusters. These points are plotted at the median redshift of each bin with 
the histogram indicating the extent of the bin. Simulations are plotted as in Figure 17. 
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The significant evolution in both P3/P0 and centroid shifts in the observations can 
be seen in Figure 18, but again this evolution is mild compared to the range of cluster 
morphologies present at any one redshift. The observed P^/Pq^ fall well within the range seen 
in our simulations at all redshifts; the main difference appears to be that in the simulations 
P3/P0 peaks at intermediate redshifts (at least for this choice of radius and luminosity cut) 
where the observed clusters have the highest P3/P0 at high redshift. The centroid shifts 
observed by Maughan et al. (2007) for clusters in the Chandra archive appear to be on 
average a bit lower than what is found in the simulations. This offset is roughly a factor of 
two, similar to the (at best) marginally significant offset seen for P2/P0 in section 3.3, and 
within the scatter in (w) nocore seen in the simulations. Since we have not used precisely the 
same derivation of centroid shifts as Maughan et al. (2007) (our definition and exclusion 
of the peak are different) nor have we reproduced all of the observational effects, we will 
not comment on the significance of this offset. However, the observations do show a similar 
upward trend of the centroid shifts with redshift with a similar slope to what is found in our 
simulations. 

In summary, on the observational side the results of Jeltema et al (2005) that clusters 
in X-ray selected samples are more disturbed at higher redshifts has been confirmed by two 
other studies using different samples and different structure measures (Maughan et al. 2007; 
Hashimoto et al. 2007), making this a fairly robust result. For a luminosity cut similar to 
that used by Jeltema et al. (2005), we do find a significant evolution in the power ratios in 
our simulations; however, we find that this evolution depends both on the cluster selection 
and radius chosen. In section 3.3, we found that at low redshift the distribution of power 
ratios in the simulations is a bit broader than observed likely due in part to our incomplete 
knowledge of gastrophysics. It is possible that this greater spread in power ratios could 
hinder detection of a mild evolution in cluster structure. We do find a significant increase 
in the average centroid shifts of clusters in our simulations with redshift independent of 
selection (R = r 500 ). The centroid shifts appear to be more robust against projection (§3.1) 
and less sensitive to deviations from hydrostatic equilibrium (§3.2.3), possibly leading to a 
more robust evolution with redshift. Evolution in a similar structure measure, the offset of 
the centroid from the surface brightness peak, was found by Kay et al. (2007) in the CLEF 
simulations, while Rahman et al. (2006) find indications of evolution in cluster ellipticity 
and multiplicity in low-redshift (z < 0.25) simulations based on the same code used here. 

We showed in section 3.2.3 that there is a strong correlation between cluster structure 
and bias in the hydrostatic mass relative to the true mass. Combined with the evolution in 
structure with redshift, though mild, this will lead to systematically larger errors in mass 
estimates at high redshift. 
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4. SUMMARY AND CONCLUSIONS 

We use Enzo, a hybrid Eulerian adaptiv e mesh refinement (AMR) /N-body code 



( INorman fc Bryan 19991 ; lO'Shea et al. 20051 ). to explore cluster morphology and its evolu- 
tion in current generation cosmological simulations. Rather than studying the impact of 
individual, idealized mergers on clusters, the large volume of these simulations allows us to 
study the development of cluster structure and its impact on observable cluster properties 
and scaling relations in a cosmological setting. We also explore the effects of projection 
on apparent cluster morphology and how well the simulations reproduce observed cluster 
structure at low redshift. We utilize and compare two commonly used measures of cluster 
structure, the power ratios and centroid shifts. Our main conclusions are summarized below. 

1) Projection: 

Based on simulations of three orthogonal projections, we find that projection along the 
line of sight leads to significant scatter in the observed structure of individual clusters. For 
example, less than 50% of clusters which appear relaxed in one projection are truely relaxed, 
and 4-10% of these "relaxed" clusters have very disturbed morphologies (for example, are 
undergoing major mergers) when viewed from another angle. 

2) Correlations with cluster properties and scaling relations: 

We do not find strong trends between cluster structure and mass, X-ray temperature, 
or Yx-i although they do correlate with centroid shift. We do, however, find a trend for 
clusters with low power ratios to have higher luminosities and a significant offset of relaxed 
clusters toward higher luminosities in the Lx — Tx and Lx — M relations. This offset is 
partially corrected by removing the cluster cores from the luminosity, but a significant offset 
still remains. When using the true cluster masses, we do not find a significant dependence 
in the Tx — M and Yx — M relations on cluster structure. 

3) Hydrostatic masses and the prospects for observable mass scaling relations: 

We find that the hydrostatic masses typically underestimate true cluster mass and that 
the error in determining mass is strongly correlated with cluster structure. For relaxed 
clusters, the hydrostatic masses are lower than the true masses by typically 12%, while for 
disturbed clusters the typical error is doubled to 24%. This correlation shows that our 
structure measures do probe deviation from hydrostatic equilibrium. We also find a weaker, 
but significant correlation of the deviation of hydrostatic mass from true mass with cluster 
mass. 

The correlation of mass underestimate with cluster structure means that when the 
hydrostatic masses are used to form the mass-observable scaling relations (Tx — Mhydro and 
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Yx—Mhydro), as would be done observationally, these relations depend significantly on cluster 
morphology. The scatter in these relations also increases. Even when using the most relaxed 
clusters to normalize the Tx — Mh y d ro and Yx — Mhydro relations, they differ significantly from 
the true relations. We can, however, use cluster structure to correct the hydrostatic masses. 
Particularly when using the power ratios, correcting the hydrostatic masses improves the 
mass estimates for 75% of the clusters, eliminates the offset between relaxed and disturbed 
clusters in the T x — M and Y x — M relations, and leads to good agreement with the true 
relations. Some difference in slope versus the true relations remains due to the dependence 
of the hydrostatic mass error on cluster mass, and the scatter in the mass-scaling relations 
is only slightly reduced. Here we do not remove substructures from our determination of 
average cluster properties, but this might allow the scatter in the mass-scaling relations to 
be reduced in high quality observations. 

The underestimate in total mass under the assumption of hydrostatic equilibrium also 
leads to an overestimate of cluster gas mass fraction. This error combines with a significant, 
but milder trend in cluster gas mass with structure. Gas mass is systematically overestimated 
for more disturbed clusters, but only by on average ~I0%. 

4) Comparison to observations: 

In this paper, we do not attempt to reproduce all of the possible instrumental and ob- 
servational effects, so we limit our comparison to observations of low-redshift cluster samples 
which have high S/N. We find that there is no statistically significant difference in either 
the average power ratios or distribution in power ratios between simulations and observa- 
tions (with the exception of a moderately significant difference in P2/P0 at small radii). Our 
simulated centroid shifts are also similar to those found by O'Hara et al. (2006), but their 
method of calculating centroid shifts is slightly different. 

Visual examination reveals that the distribution of power ratios in the simulations is 
broader than in observed clusters, including both clusters with higher power ratios and an 
extended tail of clusters with lower power ratios. Some of this difference can be accounted for 
through the addition of noise to the simulations, but it is at least partially due to incomplete 
gas physics in the simulations. In particular, our results suggest the need for a reduction in 
the peakiness of low mass cool cores but an increase in cooling in high mass clusters. The 
overall agreement between our simulations and observations indicates the general success 
of both the underlying cosmological model and the gas physics employed here, while the 
disagreements suggest that cluster structure will provide a sensitive test of "gastrophysics" 
in simulations. 

5) Evolution in cluster structure: 
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Clusters at high redshift (above z ~ 0.5), closer to the epoch of c luster formation, are 
observed to be on average more disturbed t han low redshift clusters (jjeltema et al. 2005 



Maughan et al. 20071 ; lHashimoto et al. 20071 ). In our simulations, we find a significant in- 



crease in the average centroid shifts with redshift, confirming that the evolution in cluster 
structure is also present in simulations. If we select clusters based on luminosity, similar to 
the selection of observational samples, we also find a significant evolution in the power ratios 
of our simulated clusters; however, we find that this evolution depends both on the cluster 
selection and radius chosen. 

In both observations and simulations, the evolution of cluster structure is mild compared 
to the range of cluster morphologies present at all redshifts, making its detection somewhat 
difficult and the comparison of evolution in observations and simulations sensitive to any 
differences in structure stemming from selection, gastrophysics, etc. Though mild, the evo- 
lution in cluster structure combined with the strong correlation of structure with bias in the 
hydrostatic masses will lead to an increasing error in mass estimates with redshift. 

6) Power ratios versus centroid shifts: 

Overall, we find that the power ratios and centroid shifts perform relatively similarly 
as measures of cluster dynamical state. Here we highlight a few of the areas where one 
measure performs better than another. Centroid shift appears to be in some ways a more 
robust measure. Clusters selected as relaxed or disturbed based on their centroid shifts suffer 
slightly less confusion due to projection, and the evolution in centroid shifts with redshift is 
more robust against sample selection. The power ratios, on the other hand, appear to be a 
more sensitive measure. They show a stronger correlation with deviations from hydrostatic 
equilibrium, shown in their stronger correlation to deviations in the hydrostatic masses from 
true mass. They may also be more sensitive to differences in gas physics as they show a 
slightly stronger correlation to the location of clusters on scaling relations. However, here 
we have not attempted a full comparison of these structure measures, and this comparison 
would benefit from simulations in different cosmologies and with different gas physics as well 
as further comparison to observations. 



DISCUSSION 



In terms of cosmological studies, our results are encouraging in that cluster scaling re- 
lations like Tx — M and Yx — M appear remarkably ins ensitive to cluster dynamical state . 



The biggest source of error, as noted by previous authors (IKay et al. 20041 ; iRasia et al. 2006 



Nagai et al. 2007al ). is that estimates of the total mass based on the assumption of hydro- 
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static equilibrium are systematically low. Here we also showed that the bias in hydrostatic 
mass is strongly correlated with structure. We also note that the evolution in cluster structure 
with redshift combined with the correlation of mass error with structure leads to increasing 
errors in estimates of mass and gas mass fraction with redshift. A bias exists even if only the 
most relaxed clusters are chosen to normalize the scaling relations or to measure gas mass 
fractions. In addition, our results strongly caution against the accuracy that can be derived 
from individual clusters or small samples due to the large uncertainty in cluster structure 
from projection and the very large mass errors possible for individual "relaxed" clusters. 

The goal of precision cosmology may then lead us to other means of normalizing the mass 
scaling relations, either though lensing studies (which have their own systematic biases) or, 
if deemed accurate enough, directly from simulations. For example, the fitted relationship 
in Figure 9 between the power ratios and error in the hydrostatic masses can be used to 
eliminate systematic biases with cluster dynamical state. This relationship is potentially 
testable with a study of a large sample of observed clusters of all morphologies to see if the 
normalizations of the mass scaling relations change with cluster structure as predicted in this 
paper. The offset in the simulated mass scaling relations when using the hydrostatic masses 
rather than the true masses is enough to bring them in to agreement with observations 
(see also Nagai et al. 2007b). It remains then to determine whether or not simulations are 
producing realistic estimates of the level of non-thermal pressure support in clusters. In 
the future, it may be possible to directly observe the level of bulk motions in at least small 
samples of clusters through high resolution X-ray spectroscopy or possibly the combination 
of X-ray, SZ, and lensing data (e.g. Puchwein & Bartelmann 2007). In the mean time, the 
models of cooling and feedback in current simulations are being tested in a number of ways. 
Here we showed that cluster structure can provide a sensitive diagnostic of gas physics when 
compared to observations. 

In this work, we have focused for the most part on cluster structure and its correlations 
in the simulations themselves, but in future work we will simulate mock observations of 
our clusters appropriate for the direct comparison to observed samples. Here we have used 
relatively simple indicators of cluster dynamical state based on the X-ray surface brightness 
distribution. The power ratios and centroid shift are more sophisticated than gross measures 
like ellipticity while they can still be measured fairly accurately in relatively shallow data 
(about the same depth needed for an accurate global temperature). These measures are by 
no means the only possibilities (for example, in deeper data temperature structure could 
be a sensitive measure), but we have shown here that they are effective and correlated to 
dynamical state. 

Another potentially important systematic effect of cluster structure on observational 
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studies is its effect on cluster selection in surveys. While properties like X-ray temperature 
and Yx,sz when averaged over large radii like r 500 are relatively unbiased, mergers lead 
to spikes in both luminosity and central Y$z (e.g. Ricker & Sarazin 2001; Rowley et al. 
2004; Motl et al. 2005; Poole et al. 2008). For example, the most massive cluster in 
our simulations undergoes two major mergers between z=1.5 and z=0, both accompanied by 
dramatic increases in the X-ray luminos ity near core passage when the cluster appears relaxed 



in all projections (IJeltema et al. 20081 ). In addition, cool core systems have luminosities 
which are high for their masses. The bias introduced by cluster structure on cluster selection 
in X-ray and SZ surveys is not yet well understood (nor in optical and lensing surveys which 
may be effected by projection and central concentration). 
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Table 1. Structure Measures 





R = 0.5 Mpc 




R = r 5 ao 




R = 1.0 Mpc 






Median 


Range 


Median 


Range 


Median 


Range 


P2/P0 


5.57c-6 


5.28e-7-3.20e-5 


2.56e-6 


2.70c-7-2.82e-5 


2.05e-6 


1.60e-7-2.16e-5 


P3/P0 


1.00e-7 


7.51e-9-1.03e-6 


6.65e-8 


3.58e-9-1.07c-6 


5.20e-8 


1.87e-9-1.37c-6 


P4/P0 


3.44e-8 


2.17e-9-4.28e-7 


1.92e-8 


1.01e-9-5.10e-7 


1.44e-8 


6.0e-10-4.83e-7 


(w) 


0.0359 


0.00767-0.132 


0.0192 


0.00458-0.0653 


0.00996 


0.00251-0.0304 


{^) nocore 


0.0379 


0.00852-0.138 


0.0201 


0.00505-0.0692 


0.0108 


0.00280-0.0342 



Note. — The ranges indicated in columns 3,5, and 7 show the range spanned by clusters with power ratios/centroid 
shifts greater than the lowest 10% and less than the highest 10%. The ranges are listed in this way to eliminate the 
most extreme clusters. 



Table 2. Effects of Projection 





500 kpc 


7*500 


1 Mpc 


f"500 ((w) nocore) 


% disturbed 


37.7 


30.9 


31.5 


29.6 


% disturbed appear 


55.4 


55.2 


52.7 


64.6 


disturbed in 1 proj. 










% relaxed 


16.8 


11.3 


16.8 


8.6 


% appear relaxed in 1 


10.5 


8.2 


10.6 


3.7 


proj. actually disturbed 










% appear relaxed in 1 


28.6 


38.1 


48.3 


45.0 


proj. actually relaxed 











Note. — Row 1 lists the percentage of clusters which are disturbed 
(i.e. have power ratios above our limits in at least one projection, see sec- 
tion 3.1 for definition). Row 2 lists the percentage of disturbed clusters 
which appear disturbed in any one projection. Row 3 lists the fraction 
of clusters that are relaxed (see section 3.1 for definition) in all projec- 
tions. Row 4 gives the fraction of clusters which appear relaxed in one 
projection but are disturbed in one of the other two projections, and row 
5 lists the percentage of clusters that appear relaxed in one projection 
which are relaxed in all three projections. 
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Table 3. Correlation with Cluster Properties 



Structure Measure 


M tot (r < r 500 ) 


L x (r < r 500 ) 


T x {r < r 500 ) 


Y x (r < r 500 ) 


(r < r 500 ) 


Probability 


Probability 


Probability 


Probability 


P2/P0 


0.94 


4.2e-6 


0.02 


8.7e-3 


P 3 /P 


0.69 


2.3e-5 


0.63 


0.82 


Pi/Po 


0.18 


5.9e-8 


0.70 


0.88 


yJj) nocore 


3.7e-3 


0.56 


0.030 


3.1e-3 



Note. - Probability of no correlation between each cluster property and the 
power ratios and centroid shifts based on a Spearman rank-order correlation test. 



Table 4. Mass Scaling Relations - Fixed Slope 







True Mass 
all 


all 


Hydrostatic Mass 
relaxed 


disturbed 


all 


Corrected Mass 
relaxed 


disturbed 




1 L X - E Z M 


24.959 ±0.006 


25.072 ±0.007 


25.106 ±0.014 


25.103 ±0.017 


24.962 ±0.007 


25.033 ±0.014 


24.951 ±0.017 


T x 


- E Z M 


-2.098 ±0.002 


-2.043 ±0.003 


-2.060 ±0.004 


-2.019 ±0.008 


-2.095 ±0.003 


-2.095 ±0.004 


-2.092 ±0.008 


Y x 


-E- 2/5 M 


-3.279 ±0.003 


-3.143 ±0.006 


-3.192 ±0.008 


-3.070 ±0.015 


-3.273 ±0.005 


-3.281 ±0.008 


-3.253 ±0.015 



Note. — Comparison of the normalizations of power law fits {log(y) = Alog(x) + B) to the mass scaling relations using the true masses, the hydrostatic 
masses, and the hydrostatic masses corrected based on P^/Pq. For the hydrostatic and corrected masses separate fits are also given for the relaxed and 
disturbed subsamples. The slopes have been fixed at their predicted self-similar values: 4/3 for Lx — M, 2/3 for Tx — M, and 5/3 for Yx — M. 

I 

o\ 

I 



Table 5. Mass Scaling Relations - Fit Slope 









True Mass 






Hydrostatic Mass 






Corrected Mass 








A 


B 


scatter 


A 


B 


scatter 


A 


B 


scatter 




1 L X - E Z M 


1.369 ±0.035 


24.45 ±0.50 


0.131 


1.216 ±0.036 


26.76 ±0.52 


0.152 


1.223 ±0.036 


26.56 ±0.53 


0.155 


T x 


- E Z M 


0.588 ±0.012 


-0.96 ±0.18 


0.078 


0.527 ±0.014 


-0.04 ±0.21 


0.112 


0.529 ±0.014 


-0.11 ±0.21 


0.110 


Y x 


- E Z 2/5 M 


1.601 ±0.022 


-2.39 ±0.32 


0.045 


1.435 ±0.032 


0.12 ±0.46 


0.085 


1.442 ±0.032 


-0.09 ±0.46 


0.082 



Note. — Comparison of the power law fits (log(y) — Alog(x) + B) to the mass scaling relations using the true masses, the hydrostatic masses, and the 
hydrostatic masses corrected based on Ps/Pq- These fits include all clusters, and the slope is free to vary. Scatter is measured as Ei(^°.9( x i) + B/A — 
log{ yi )/Af/NY/\ 

I 

I 
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Table 6. Comparison to Low Redshift Observations 





Median Sims. 


R = 500 kpc 
Median Obs. 


Rank- Sum 


KS 


Median Sims. 


R = 1 Mpc 
Median Obs. 


Rank- Sum 


KS 


P2/P0 


4.4e-6 


2.1e-6 


0.02 


0.002 


2.6e-6 


1.0e-6 


0.10 


0.07 


P3/P0 


4.3e-8 


3.4e-8 


0.40 


0.63 


3.0e-8 


5.3e-8 


0.34 


0.14 


P4/P0 


2.1e-8 


l.le-8 


0.17 


0.04 


1.3e-8 


9.2e-9 


0.33 


0.12 



Note. — Comparison of the power ratios of our simulated clusters to the ROSAT observed sample from Buote & 
Tsai (1996) at aperture radii of 0.5 Mpc and 1 Mpc. Columns 2 and 6 gives the median power ratios of the simulations 
compared to the medians of the observations in columns 3 and 7. Columns 4 and 8 give the probability from a Wilcoxon 
rank-sum test that the mean power ratios of the observed and simulated samples are the same. Columns 5 and 9 give 
the probability from a Kolmogorov-Smirnov test that the distributions of power ratios are the same. 



Table 7. Evolution 







R = r 5Q0 






R = 0.5 Mpc 






Low vs. Mid 


Mid vs. High 


Low vs. High 


Low vs. Mid 


Mid vs. High 


Low vs. High 


Mass cut 














P2/P0 


0.20* 


0.061 


0.18 


0.31* 


0.36* 


0.19* 


P3/P0 


0.40* 


0.14 


0.17 


0.18 


0.28 


0.10 


P4/P0 


0.15* 


0.077 


0.24 


0.38* 


0.26 


0.30 


{' w )nocore 


0.44 


0.011 


0.0053 








Lx cut 














P2/P0 


0.09 


0.16 


0.0092 


0.10 


0.015* 


0.18* 


P3/P0 


0.0043 


0.44 


0.0023 


0.0093 


0.13* 


0.13 


P4/P0 


0.048 


0.23 


0.0058 


0.058 


0.12* 


0.28 


(lv)nocore 


0.0036 


0.23* 


0.023 









Note. — Probability based on a rank-sum test that the means of each structure measure are the same 
between the specified redshift bins. Three redshift bins are considered: low rcdshifts (z < 0.25), intermediate 
redshifts (0.25 < z < 0.65), and high redshifts (0.65 < z < 1.5). Probabilities are tabulated for clusters 
selected based both on a mass cut of M > 2 x 10 14 M Q and a luminosity cut of Lx > 2 x 10 44 ergs s _1 
(*) Denotes that the mean of the structure measure decreases between redshift bins rather than increasing. 



